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Abstract 

We seek analytic formulae for the pricing of non-payment insurance on a vanilla interest rate swap. 

In the event of a default bv the swap counterparty, the net PV of all future cash flows, floored at a given 
strike level (typically zero), is required to be paid by the insurer to the insured party at the time of 
default. We suppose that the domestic rates and the counterparty default intensity are stochastic and 
potentially correlated. The credit intensity model is taken to be Black-Karasinski and the interest rate 
model to be beta-blend, a family of models which encompasses within its scope both the Black-Karasinski 
and the Hull-White cases. A Green’s function solution for the governing PDE is found as a perturbation 
expansion valid in the limit of low interest rates and credit default intensities. This is used to calculate 
to first order accuracy the fair price of insuring payer or receiver swaps, including the impact of capping 
the insurance protection amount. 

1 Introduction 

We seek analytic formulae for the pricing of non-payment insurance on a vanilla interest rate swap. We 
suppose in the first instance that the insurance buyer receives floating rate payments from the swap 
counterparty and pays fixed. In the event of a default by the swap counterparty (floating rate payer) 
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before some specified time T, the net PV of all future cash flows, floored at a given strike level (typically 
zero), is required to be paid by the insurer to the insured party at the time of default. The approach taken in 
the calculation is similar to that of Turfus (2017a) where the cost of insuring an equity option was considered. 
Use is made below of the methodology set out in the work of Horvath et al. (2017) and Turfus (2017b). 

We start in section 2 by defining the underlying processes in terms of suitably chosen auxiliary variables 
and the no-arbitrage conditions they are required to satisfy; we infer therefrom the governing PDE for 
contingent claims. Given the intractability of this PDE, we define in section 3 an asymptotic scaling of the 
underlying variables and associated functions based on the assumed smallness of the two short rates: interest 
rate and credit default intensity. This gives rise to a tractable leading order PDE with the “intractable” parts 
isolated as a perturbation: we can then seek solutions as perturbation series. Rather than looking to obtain 
particular solutions directly, a Green’s function for the full PDE is sought as a perturbation expansion. 
The constituent cash flows associated with the interest rate swap are valued conditional on a default at 
time t < T in section 4 in terms of the asymptotically scaled variables. The sum of these contributions, if 
positive, is the exposure payment which must be made by the insurer. The further calculation required to 
derive the fair value of insurance for each of the constituent conditional cash flows is presented in Appendix 
A. The contributions are summed and the aggregate price of insuring all cash flows associated with payer 
and receiver swaps is presented in section 5, including specialised results valid in the (singular) Hull-White 
limit. Consideration is also given to assessing the impact of a cap on the insurance payment amount. Finally 
the predictions of the asymptotic model are compared with results obtained from a Monte Carlo simulation 
in section 6. 


2 Stochastic Modelling 

2.1 Definition of Underlying Processes 

Our modelling approach is to represent the interest rate and the credit default intensity as short rate diffusion 
processes. Specifically, we suppose the interest rate process ry to be governed for t > 0 by a beta-blend short 
rate model as set out in Turfus (2016), Horvath et al. (2017) and Turfus (2017b). This encompasses both 
the Hull & White (1990) and Black & Karasinski (1991) short rate models. We suppose the credit default 
intensity process X t to be governed by a Black-Karasinski (lognormal) short rate model. 

We shall find it convenient to work with auxiliary variables ay and y t satisfying the following Ornstein- 
Uhlenbeck processes: 

dx t = — a r x t dt + <J r (t) dWt, 
dy t = -a x y t dt + a x (t) dWf, 

where dW^ and dW? are Brownian motions under the usual equivalent martingale measure with 

COTT(W t 1 ,W?) = PrX 

assumed and Xq = yo — 0. These auxiliary variables are related to the interest rate ry and the credit default 
intensity X t , respectively, by 

(1 - P)r t + pr{t) = ( r(t ) + (1 - /3 )r*(t))£ , (2.3) 

X t = (X(t) + X*(t))£(y t ), (2.4) 

with /3 £ [0, 1) assumed. Here r(t ) is the instantaneous forward interest rate and A (t) the associated credit 
spread (see (2.9) below), with a r (t) and u x (t) their respective volatilities and a r and a x their mean reversion 
rates. Here £{X t ) := exp (Xj — i[A] t ) is a stochastic exponential with [A] t the quadratic variation of a 
process X t under the requisite measure. We note here specifically that \x\ t = I r (t ) where we define 

I r (t) := f e _2ar ( t_ “Vr(u) du. (2.5) 

Jo 


( 2 . 1 ) 

( 2 . 2 ) 
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The required form of the configurable functions r*(t) and A*(f) is determined by calibration of the model 
to satisfy the no-arbitrage conditions set out below. We further assume that Xq = yo = 0, with t = 0 the 
“as of” date for which the model is calibrated. As can be seen, the beta-blend model for the interest rate 
represents a hybrid between the Hull-White model ((3 f 1) and the Black-Karasinski model (/? = 0). 

The no-arbitrage conditions 

The formal 110 -arbitrage constraints which determine the functions r*(t) and A *(t) are as follows. First, by 
considering a risk-free cash flow at time t, we deduce 

E r e -/o r * ds l =D{0,t), (2.6) 

for 0 < t < T m , where T m is the longest maturity date for which the model is calibrated, and 

D{t 1 ,t 2 )=e~^ Hs)ds (2.7) 

is the fi-forward price of the ^-maturity zero coupon bond. Likewise, considering a risky cash flow, we have 

E e -/o( r »+ A <>) ds =B(0,t), (2.8) 

where _ 

B(ii,t 2 ) = e ~ f ti 2(r(s)+A(s))ds (2.9) 

is the fi-forward price of a risky (defaultable) cash flow. We shall assume the bond prices can be ascertained 
at the initial time t = 0 from the market, whence we can view (2.7) and (2.9) as implicitly defining r(t) and 

A (t). 

2.2 Derivation of Governing PDE 

We consider the problem of pricing a security with maturity T which, in the event of a counterparty default 
at some stopping time r € (0, T], pays a (domestic currency) cash amount which depends on f := x T (but 
notably not on y T ), denoting this payoff by P(£, r). We further introduce the convenient shorthand notation 
that, for a process X t and deterministic function /(.), 

£x(f(t)X t )) :=£(f(t)X t ))\ Xt=x , 

in terms of which we can re-write (2.3) and (2.4) as rt = r(xt,t) and X t = X(y t ,t), where 

r(x,t) := ^(r(f) + (1 - /3)r*(t))£ & ^ ^ - /3r(f)^ , (2.10) 

A (y,t) := (X(t) + X*{t))£y(y t ), (2.11) 

Writing the price of the security at time t £ [0, r A T] as f t = f(xt, jjt,t ), we can infer by application of 
the Feynman-Kac theorem to (2.1) and (2.2) in the standard manner that the function f(x,y,t) satisfies the 
following backward diffusion equation: 

(J^ + t - r(t) - A (t)j /(:r, y, t ) = -A (y, t)P(x, t), (2.12) 

subject to f(x,y,T) = 0, where 

£ := - a '*§& - + l + +^(‘)|i) - «*.<) - ? (*))- 

- (A(h, 8) - A(e))- (2.13) 

In the absence of closed form solutions to (2.12) and guided by the work of Hagan et al. (2015), Horvath et 
al. (2017) and Turfus (2017b), we propose a perturbation expansion approach as follows. 
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3 Asymptotic Modelling 


For both short rate models we formally apply a “low rates” assumption. To this end we define small 
parameters 


and 0(1) functions 


and 


1 


Oi r Tm J o 

i r 1 


■= 


a\T r 


"(t) dt 


\(t) dt 


m J 0 


r*(t);= e r l r*(t) 

H*) = 


(3-1) 

(3.2) 

(3.3) 

(3.4) 


Note, we make no formal requirement that either oy(t) or cr\(t) be small (scaled, say, in terms of their 
respective mean reversion rates). But, if they are, this will in general improve the accuracy of the perturbation 
expansions deduced. We further define a new scaled variable xt and associated scaled volatility cr x (t) by 


x t := e r ^x t e art , 
0x(t) = e~ p a r {i)e art . 


Likewise we define 

Vt ■■= Vte a ^, 

v y {t) ■= cr\{t)e axt . 


(3.5) 


(3.6) 


Here the exponential time scaling is to facilitate removal of the mean reverting drift terms in (2.13). We 
further define new functional forms /(•) and P(- ) by: 


f(xt,Vt,t ) := e* 1 f(x t ,y t ,t), 
P(x t ,t) := P(x t ,t), 


where x t and y t are related to Xt and y t by (3.5) and (3.6), respectively for t £ [0, r A T\. In this notation, 
(2.1) and (2.2) can be re-expressed as 

dx t = cr x (t)dWf, (3.7) 

dy t = cr y (t) dW? (3.8) 

and (2.12) as 

(J^ +C- f(t) - A (t) - e r h( x, t ) - e\g(y, f) j f(x, y, t) = -(A (t) + g{y , t))P{ x, t) (3.9) 

where 

£ N := 2 (^^dx? + ^ PrX<Jx ^ ay ^Qxdy +(Ty ^dip) (3.10) 
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is a pure diffusion and 


h(x, t) := h(x, t, t), (3.11) 

g(y,t) ■■= g(y,t,t), ( 3 . 12 ) 


with 


h(x, t, ti) 
Fp(t) 

g(y,t,h) 


+ f*(ti))£ x (Fp(ti)x t ) - r(t i), h > t, 
(1 - pY-Pe- 0 "- 1 
I r{t)\P 

(A(ti) + X*(t 1 ))£ y (e~ axtl y t ) - A(ti), h>t. 


(3.13) 

(3-14) 

(3.15) 


We seek a Green’s function solution for (3.9) as a joint power series in e r and e\. The reader is referred 
to Turfus (2017b) for details of the derivation up to second order accuracy. We note here the important 
observation there made that satisfaction of the no-arbitrage conditions requires that r*(t) ~ e r r\{t) and 
A *{t) ~ e r X{flit) + e AA^i (t) with errors 0(e ji) and 0(e % + e\), respectively, where 


with 


r i / r(u) (exp(Fp(u)Fp(t)I x (0,u)) - l)du, 

Aq i(t) = A (t) j A (u) ^exp ^e _aA ^ t+1i ^/j / (0, u)^ — l^j du. 

Xi,o(t) = r(t) [ X(u) (exp (e~ axu Fp(t)I xy (0,u)) - l) du 

+ X(t) [ f(u) (exp (e _aAt F /3 (w)7 X j / (0,u)) - l) du, 

Jo 

rt 2 

Ix{ti,t 2 ) ■■= / al(u)du 
Jt i 
rt2 

Iy(ti,t 2 ) := / &y(u) du. 

Jti 

CT x {u)(Jy{u) du. 


Other details we shall report as needed. 


(3.16) 

(3.17) 

(3.18) 

(3.19) 

(3.20) 
(3-21) 


4 Exposure Modelling 

As stated in §2.2 above, the insurance payoff is determined by the expected value of all future cash flows 
associated with the derivative insured at the time of a counterparty default, floored at zero; in other words 
the positive exposure. These flows are of two types: 

1. fixed coupon flows, and 

2. Libor flows. 

To obtain an expression representing P(x,t) in (3.9), we need to evaluate these cash flows at arbitrary times 
t £ (0, T), independently of whether the counterparty defaults or not. 
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4.1 Conditional Pricing of Component Flows 

Fixed Coupon Flow 

We consider the PV as of time t < T of a fixed flow of one unit of domestic currency at time T, denoting 
this as Vp(xt,t). Turfus (2017b) finds that 

Vp(x,t) ~ D(t,T) (1 + e r Fp 0 {x,t) + e 2 r Fp fi (x,t)) (4.1) 


with 0(e)!) errors, where 

Fp 0 (x,t) := - j r(u)(£ x (Fp(u)x t ) - 1) du 
F£ 0 (x,t) := + J £ x (Fp(v)x t ) 


f{v) / r(u)£ x (Fp(u)x t ) (exp(F /3 (n)F /3 (n)4(0,M)) - 1) du - r*(u) ) dv. 


We shall make use here only of the first order expression, writing 

Vp(x,t) = D(t,ti ) (l + e r F 1 ‘( 0 (x,t)) + 0(e)!), t < U. 


(4.2) 


Libor Flow 

Let us next denote the PV as of time t < U of the Libor flow fixing at ti_i and paying at time ti by Vpf 1 (x t ,t) 
per unit notional. This can be expressed as 


Vp 1 (x, t) - Vp (x, t), t < ti_i 

Vp(x,t) (Vp(x,ti- 1) _1 - 1) , U -1 <t<U, 


(4.3) 


the latter expression being required to calculate the impact of a future Libor payment which has already 
fixed at the time of counterparty default. Noting that U) -1 — 1 is in fact first order in e r , the value 

of the Libor payment can be rewritten as 




1 ) (l + e r Fp 0 ^aqf)) - D(t, U) (l + e r F 1 t ( 0 (a:, t)) +0(e 


t < ti-i 


DiU-ut )- 1 1+erFl ’°( x,t ) 


1+^,;;;^) -£>(*, *<)(i + ^(x, *))+ o(e% < t<u. 


(4.4) 


4.2 Aggregating the Payoff Flows 

Next we compute the positive exposure function P(x,t) appearing in (3.9) above using the expressions for 
the component flows derived in the previous subsection. We consider for simplicity that the coupon flows 
occur at the same dates and the day count convention used is the same on both legs. 1 Suppose that the 
payment periods are fj], for i = 1 ,..., n and the day count fractions A*. Consider a payer swap where 
the pay leg has a fixed coupon amount c and the other leg receives Libor. 

In terms of the above notation the value of the positive exposure at default time r < T is then given by 

P(x,t) = max | Y. (Vl\x,t) - cA iVp(x,T)^ l 7 -<min{t i ,T}, oj (4.5) 

1 I(. is not difficult to modify the formulae derived to take account of a difference in day count convention between the two 
legs. 
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per unit notional. It remains to solve (3.9) for f(x,y,t) subject to (4.5) and a homogeneous final condition 

f(x,y,T) = 0, (4.6) 

to the desired order of accuracy in the parameters e r and e\. In practice we will restrict our attention to 
first order in each, to which end we will need only the leading order representation of the Green’s function 
derived in Turfus (2017b), in combination with first order versions of the expressions derived in section 4.1 
above. This calculation is performed in Appendix A below, the main conclusions of which are set out in the 
following section. 


4.3 Aggregated Payoff for Hull-White Case 

For future use, we express the exposure expression in (4.5) in original unsealed variables for the particular 
case of a Hull-White interest rate model. We obtain straightforwardly, in self-evident notation, 

n 

Exposure = ^ (v^ (x, r) - cAiVp (x, r)) , (4.7) 

i=1 

where, using the first order approximation for the fixed flow values as of time t < ti, we have 

Vp (x, t) ~ D(t,ti ) (1 - xB*(ti - t )), (4.8) 

with errors = O(e^), where 

1 _ e ~a r v 

B*(v):= -. (4.9) 

Qlf 

Likewise we have to the same level of accuracy 




D(t, U- 1 ) (1 - xB*(U-i - t)) - D(t, U) (1 - xB*(U - t)), t < U-i 
IHJi -T- 1)- 1 ~ I 1 - xB*(U - t)), U-1 < t < U. 


(4.10) 


5 Main Results 

PV Insurance for Payer Swap 

We conclude from the asymptotic analysis that, with relative error = 0(e r + e\), the cost of insurance 
purchased at t = 0 on a payer swap is given in original notation by 

n 

insurance ~ ~ C fp A >) ’ (S' 1 ) 

i= 1 


where 


pwAT 

fp:= / \(v)B(0,v)D(v,w)N(-d 1 (£*(v),v))dv 

Jo 

rtiAT 


(5.2) 


cti r(u) 


Sl : = ( D (ti- iA) _1 - 1) fp +D(t i - 1 ,t i )- 1 J A (v)B(0,v)D(v,ti) J —_ ^ 

(exp ti-i)I r \(v A ti-0^ N (-d 2 (£*(«)> u, v, v A U-i)j 

— N (— di(£*(v), v)j )dudv , (5.3) 
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with 


I r x{t) ~ Pr\ [ e (“r+<**)(t U ^a r (u)ax(u)du, 

Jo 


di(x,v) := 


x - I r x(v) 

\/lr{v) 


d, 2 {x, u, v, w) := d\(x, v) — 


Fp{u, 2w — v)I r (w) 

A(*') 


(1 _ p) e ~ ar< ' u ~ v '> 
C(v) e~ arV x* (v), 


(5.4) 

(5.5) 

(5.6) 

(5.7) 

(5.8) 


where x*(v) is given by (A.2). Alternatively, in the Hull-white case, f*(v) is given to leading order as the 
zero of (4.7) with r = v, subject to (4.8) and (4.10). 

Note that in the Hull-White case (/3 = 1), we must take the limit as /3 —► 1 in the inner integral in (5.3). 
Applying l’Hopital’s theorem and computing explicitly the inner integration w.r.t. u, we obtain: 


/i/ _• y 1 — 1 ) fp + ti) 


rUAT 


-l 


\(v)B(0,v)D(v,ti)B*(ti - U- 1 ) 


7(^-i) v)I r x{v A (v), v)) + e a ’' |u u -^I r (v A t i - 1 ) N ^ 1 dv, (5.9) 

V V J rW / 


with 


7 (u,v) := 


' e -<xx(v-u) t „< W) 
e -a r ( U - V ) u>v . 


(5.10) 


This completes our leading order calculation. We note in parting that the sum of fp in (5.1) can be computed 
more efficiently in practice by refactoring as 


*=i i=i 


rtjAT 


\(v)B(0, v)D(v, tj)N(—di(f* (v), v)) dv, 


(5.11) 


i= 1 


tj — iAT 


where we define to = 0. The first term in (5.9) can be handled identically. As the integrand is now 
independent of i, it effectively only needs to be integrated once, from 0 to T. 


PV Insurance for Receiver Swap 

The corresponding expression for the receiver swap is easily derived. (5.1) is replaced straightforwardly by 




insurance 


n 


i=1 


(5.12) 


with the caveat that —dj(-) is replaced by dj(-) in (5.2) and (5.3). In the Hull-White case we have 


f ( L ~ 1 - l) fp + D(ti-i,ti) 


-1 


^min{^i,T} 


A (v)B(0,v)D{v,t i )B*(t i -t i - 1 ) 


liti-i, v)I r \(v A ti-i)N(di(£*(v), v)) — e arlv u -^I r (v A U-i) N ^ ) dv. (5.13) 

V \Aw) / 
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PV Insurance with Capped Payment 


We consider finally how to handle the situation where the insurance payment is capped at a given level, K 
say. We deal with the case of a payer swap but the calculation is exactly analogous for a receiver swap. Also, 
although we assume that K is constant for notational convenience, it is a straightforward matter to modify 
(5.14) and (5.15) below to take account of a time-dependent cap. First let us consider the modified form of 
the payoff if the critical boundary were at K rather than 0, viz. 


Pk(x, t) = max 


( V L l) ( x > T ) - c Wf Ou r )) 1 r<t i - K, 0 


(5.14) 


Further, let us denote the critical value of x at which this payoff comes into the money by x* k (t) and the 
corresponding value of x by &(r). Let the corresponding PV of fixed and floating payments paid conditional 
on this critical value being exceeded upon counterparty default be written as fp K and f L K , respectively. 
It remains to value the additional effective payment of K at default in this version of the contract. The 
calculation is almost identical to that for a fixed cash flow at T except in that the cash flow is instead brought 
forward to r. We obtain the following leading order estimate: 2 


fx = K / X(v)B(0,v)N(-d 1 (C K (v),v)) dv 


(5.15) 


and conclude 


V- K 

v insurance 


£(/& P//.1/A) / 


T 

F.K i 


(5.16) 


Finally we observe that the required payoff for an insurance payment capped at K is P(x,t) — Px{x,t), 
whence, by the principle of portfolio replication, the PV of the capped insurance will be 


Vi 


ins_capped 


= K 


insurance 


— V: 


K 

insurance * 


(5.17) 


6 Comparison with Monte Carlo Results 

Comparisons were made between prices based on our formulae and those generated by a Monte Carlo 
representation of the same model for the following base situation. Parameters were varied in the manner 
specified to see their impact on pricing and on the accuracy of our approximation. All prices below are 
reported in M EUR. 


Product details 

Swap maturity 

Swap coupons 

Notional 

Protection period 
Insurance coupon rate 

lOy 

Receiving quarterly 2.5% fixed; paying 6M LIBOR + 40bp spread 

100M EUR 

6 y 

Paying semi-annually 4% 

Market details 

Interest rate 

Short rate rol (normal) 
Market CDS rate 

Credit vol (lognormal) 
Correlation 

Swap rates increase from -24bp to 81bp during the lifetime of the swap 

cr r (t) increases from 0.17% to 0.81% during the insured period with a r = 25% 

Flat 3.5%:. with recovery rate 55% 

a\(t) = 70% with a\ = 30% 

p r \ = 30%. 


2 If we consider K to be 0(e r ), the level of accuracy of this expression will be the same as those previously computed. But 
even if in practice K is chosen to be relatively large, this will usually not be too much of a problem since the payoff associated 
with the cap will then typically be significantly out of the money and so constitute only a small correction to the uncapped 
value. 
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The impact of interest rate volatility in the absence of correlation is illustrated in Fig. 1. Here we expect 
no impact of credit volatility. A relative shift of the entire volatility curve is used: -100% corresponds to 
no volatility, while +100% represents a doubling of the volatility levels. As can be seen, the asymptotic 
model predictions follow the Monte Carlo results closely. However it is notable that there remains a residual 
difference between the models even when the volatility is zero, at which point the asymptotic model becomes 
exact. This difference is a consequence of an approximation used in the Monte Carlo model, which does not 
model default events continuously but at discrete weekly intervals. This approximation results in a small 
underestimate of the value of protection. We conclude that our asymptotic formula is calculating the PV of 
the insurance correct to within a basis point of notional. 
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Figure 1: PV dependence of interest rate swap insurance on interest rate volatility level. 

The impact of correlation on prices is illustrated in Fig. 2. The small discrepancy of around 1 bp of 
notional at zero correlation is again evident. The agreement for non-zero correlation is generally good, at all 
times within 6 bp of notional even in the most unfavourable case with a correlation of -100%. But there is 
a clear trend for the asymptotic formula to predict a greater impact of correlation than is suggested by the 
Monte Carlo simulations. This difference may be due to higher order terms which have not been calculated, 
specifically at 0(e r e\). More detailed calculation is required to resolve the issue. 
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Figure 2: PV dependence of interest rate swap insurance on rates-credit correlation level. 

Finally the impact of credit volatility is examined in the presence of a 30% correlation (since there is no 
impact in the absence of correlation). Results are illustrated in Fig. 3. The agreement is seen to be excellent 
for credit volatilities up to a rather large value of 80%. The slight discrepancy between the gradients of the 
two graphs is attributable to the same cause as in the previous case considered. 
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Figure 3: PV dependence of interest rate swap insurance on credit volatility level. 


A Asymptotic Pricing 

Insurance Pricing 

To solve (3.9) for f(x,y,t ) subject to (4.5), we must first determine the region in (x,t) space within which 
the exposure is positive. Since the exposure is clearly a monotonically increasing function of x. which is 
furthermore negative in the limit as x — > —oo but positive for large enough x, we infer, for fixed r, the 
existence of a critical value of x above which the exposure will be positive. In other words, there exists x*(t) 
such that 

x > x*(t) <t=> P(x,t) > 0. 

The positive exposure function in (4.5) can on this basis be rewritten as 

n 

P{x,t) = (v^\x,t) -cAiVp(x,r)j 1 t <u,x>x*(t)- A 1 ) 

i=i 

If we denote by Pq(x,t) the approximation to the r.h.s. of (4.5) obtained by substituting in the first order 
approximations (4.2) and (4.4), we can calculate x*{t) to leading order for a given r as 

x*{t) ~ inf {a; | Pq{x,t) > 0} (A.2) 

with errors = 0(e r ). We will look therefore to calculate /(•) correct to 0(e r c\), making use of the results 
in §4.1. The linearity of the positive exposure function and the fact that f(x,y,t) satisfies the homogeneous 
final condition (4.6) then allow us to obtain the fair price of insurance by solving (3.9) term by term with 
the help of our asymptotic Green’s function. In fact we shall only be interested in the value of /(•) at t = 0 
(with also x = y = 0), so we shall make this simplification from the outset. 
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Difference 
































Thus we solve (3.9) asymptotically for t = 0 with P(x,t) taken in turn to be each of Vp\x,r) and 
Vp(x,T ), contingent on x > x*(t), denoting the results by and fp , respectively. The aggregate cost of 
insurance is then given by 

n 

/( 0 , 0 , 0 )~^(/«-cA i /^) (A.3) 

i—1 


Fixed Flow 


Starting with the fixed leg and making use of the expression (4.2), we see that calculation of fp to leading 
order requires use only of the leading order Green’s function 


d 2 

G o,o0 v,y,t;t,T),v) =B(t,v)-^ K: N 2 (t-x,ri-y,R(t,v)), t< 


d^dy 


(A-4) 


where N$(x, y, z\ R(t , v)) is a bivariate Gaussian probability distribution function with mean 0 and covariance 
matrix 

I xy (t,v) 

Ixy{t,v) Iy(t,v) 

Setting t = 0 and carrying out the integration over £ and 77 , we obtain 


R(t, v) := 


(A.5) 


pOO pOO 


n 


B(0,v)D(v,ti)\(v) 


'X* ( V) J — OO 


£ v (e~ axv y v 


d 2 


d^dy 


y; R( 0, v)) dyd( dv 


B( 0, v)D(v, ti)\(v)N(—di(x*(v), v, v)) dv, 


(A. 6 ) 


with relative error = 0(e r + e\), where 


di(x,v) := 


x-e axv I X y(0,v) 

\Z J x( 0,v) 


(A-7) 


Libor Flow 


Considering next fp\ we see that from the form of (4.4) that this has two contributions, one implicitly 
0(e r ) and the other explicitly so. Let us therefore write 


f(«) 

JL 


f{i) 
J LI 


f(») 

rJ L 2 


(A. 8 ) 


in obvious notation. The first term is, for a given term structure of interest rates, equivalent to a fixed flow 
at ti given by the forward Libor rate, whence we deduce 

fLi~(D(ti-uti)- 1 -!)^ (A.9) 

with relative error = 0(e r + e\). For the other term, we have 


f(i) 
J LI 


poo poo 


D(ti— 1 , t) 1 / B(0,v)D(v,ti)X(v) 


& 2 


f {£ x (Fp(u)x vMi -i) - l) r(u) dud^dydv 


1 / B(0,v)D(v,ti)X(v) / r(u) 

do Jti -1 

(exp [Fy{u)e~ axV I xy {Q, v A tj-i)) N (—d 2 {x*(v),u, v, v A — N (—di(x*(v), v))) du dv, (A. 10) 


where 


d^{x, u, v, w) 


di(x,v) 


Fp(u)I x (Q,w) 

\/lx{ 0,u) 


(A.ll) 
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